Unstable superheavy relic particles as a source of neutrinos responsible for the 

ultrahigh-energy cosmic rays 
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Decays of superheavy relic particles may produce extremely energetic neutrinos. Their annihilations 
on the relic neutrinos can be the origin of the cosmic rays with energies beyond the Greisen-Zatsepin- 
Kuzmin cutoff. The red shift acts as a cosmological filter selecting the sources at some particular 
value z e ± Sz, for which the present neutrino energy is close to the Z pole of the annihilation 
cross section. We predict no directional correlation of the ultrahigh-energy cosmic rays with the 
galactic halo. At the same time, there can be some directional correlations in the data, reflecting 
the distribution of matter at red shift z = z e ± Sz. Both of these features are manifest in the existing 
data. Our scenario is consistent with the neutrino mass reported by Super-Kamiokande and requires 
no lepton asymmetry or clustering of the background neutrinos. 
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Cosmic rays [jj] with energies beyond the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff ^ present a challenging 
outstanding puzzle in astroparticle physics and cosmol- 
ogy. One of the most appealing and economical explana- 
tions, proposed by T. Weiler |Jj, is based on the obser- 
vation that the GZK cutoff is absent for neutrinos HQ . 
High-energy neutrinos can annihilate near to Earth on 
the background relic neutrinos via the Z resonance || 
producing the ultrahigh-energy cosmic rays (UHECR). 
A very appealing feature of this scenario is that the req- 
uisite energy scale appears naturally p] if the relic neu- 



trinos have mass m v = m c 



O.Q7±8;S eV, 



as inferred 



from the Super-Kamiokande (SK) data (assuming no de- 
generacy in the neutrino mass spectrum) . Such relic neu- 
trinos can efficiently annihilate a neutrino with energy 
E v = M 2 /2m SJf ~ 10 23 eV, corresponding to a Z res- 
onance. This energy happens to be just a few orders of 
magnitude above the GZK cutoff, which is precisely what 
is needed to explain the GZK puzzle. 

The only problematic element in this scenario |||| so 
far is the source of the ultrahigh-energy neutrinos. In 
this paper we show that heavy unstable relic particles can 
produce the requisite neutrinos. The spectrum of these 
neutrinos is different from that of known astrophysical 
neutrino sources (it grows rapidly with energy, up to a 
sharp cutoff), so that the upper bounds of Refs. || do 
not apply. The flux of ultrahigh-energy neutrinos may be 
high enough to explain UHECR regardless of the lepton 
asymmetry. 

Heavy relic particles decaying into hadrons and pho- 
tons have been considered as sources of UHECR |Q,pl?,|| . 
However, the lack of directional correlation of the ob- 
served events with the galactic halo and the observed 
clustering of events pi in the directions unrelated to 
the distribution of the galactic dark matter may stymie 



this explanation. Indeed, if the decaying particles are 
the origin of the ultrahigh-energy protons and photons, 
one expects the directions of UHECR to point to some 
nearby sources. The present data are too sparse to reach 
a definite conclusion [jll| . 

If, however, the superheavy relic particles decay into 
neutrinos, the directional correlation with the local distri- 
bution of matter is not expected in general. The only sig- 
nificant contribution to the vv annihilation (vv — > Z — > 
P7-..) on background neutrinos with mass m v comes from 
energies near the Z pole, so the high-energy neutrino 
must have energy close to E Tes — M|/2m„ to produce 
an appreciable amount of cosmic rays. Let us suppose 
that the ultrahigh-energy neutrinos come from a two- 
body decay of some superheavy particle X into a neu- 
trino and some other light particle. If the mass of the 
relic particle m x is not fine-tuned to be exactly 2E rcs , 
the local distribution of A-particles has no effect on the 
highest-energy cosmic rays. If m x > 2E rcs , the neutri- 
nos produced at red sift z e = (m x /2E rcs — 1) annihilate 
very efficiently close to Earth because their present en- 
ergy corresponds to the Z pole. The red shift acts as a 
cosmological filter selecting the sources at some particu- 
lar value z e . The directions of the UHECR should, there- 
fore, point to sources at z = z e , not to those in the local 
group. One may still observe clustering of events conge- 
nial to the distribution of luminous matter at z ~ z e , but 
not at z w 0. 

Let us consider some superheavy relic particle X that 
can presently contribute to the cold dark matter p|,p^[. 
Such particles could, for example, be produced non- 
thermally at preheating M. Let us also assume that it 
has a lifetime r x that is much greater than the age of 
the universe to and that it decays with emission of at 
least one neutrino with branching ratio B v . We assume 
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that hadronic and other visible decay modes have a neg- 
ligible branching ratio. For simplicity, we will consider a 
two-body decay with emission of a neutrino with energy 
E u , e — m x/2- Generalizations to many-body processes 
and other modes of decay are straightforward. 

The spectrum of ultrahigh-energy neutrinos from de- 
caying unstable heavy particles is given by p"3| . |l~i| 

where we have assumed t x 3> to- Here </> 7! o = ri 7j o/47r 
10 12 cm _2 s _1 sr _1 is the flux of the cosmic microwave 
background photons, and Y x is the ratio of the X num- 
ber density to that of photons, Y xo — n xo /n 7 .o- We 
assume that the universe is matter-dominated from the 
time of emission until present. As usual, index refers 
to the present time, while subscript e denotes the time 
of emission. 

We will now compute the flux of UHECR from the 
annihilation of ultrahigh-energy neutrinos on the back- 
ground neutrinos with mass m v . The cross section for 
the process vv — ► Z* — ► hadrons has a sharp maxi- 
mum at the Z pole: a oc 1/ {(s - A/ 2 ) 2 + M 2 T 2 }, where 
s = 2E u m v , and T z sa 2.49 GeV is the Z width. We are 
interested in the neutrino annihilations within the atten- 
uation length of protons and photons, which is about 
50 Mpc from Earth. Within this radius one can neglect 
the red shift. Only those neutrinos that have energies 
E vfi = E rcs ±SE = (M*/2m v )[l ± T z /M z ], can con- 
tribute to the flux of UHECR. The flux A<f> of these 
neutrinos can be obtained by integrating d(f>(E)/dE in 
equation ([[]) from E IeB — SE to E IeB + SE: 

The value of n x = fl x p c /m x is bounded from above 
by the requirement that Sl x h 2 < 0.25. X particles can 
be the cold dark matter §,[l| if fl x h 2 « 0.1 - 0.25. 

The probability P for each of these neutrinos to in- 
teract inside the GZK sphere, that is within 50 Mpc 
from Earth, depends on the lepton asymmetry param- 
eter rj = ?ii/,relic/^7,o- We treat rj as a parameter whose 
value is constrained by nucleosynthesis, as well as large- 
scale structure and cosmic microwave background (CMB) 
radiation anisotropy measurements. 

If the universe has zero lepton asymmetry, 77 = 0.14. 
However, depending on the cosmological parameters, this 
may not be the value that gives the best agreement with 
observations. If the cosmological constant is small or 
zero, 51a < 0.1, the value 77 w 4 is favored by the 
combination of the large scale structure and the CMB 
data For Q A w 0.5, the range < 77 < 2 is al- 

lowed [ft6LFor larger values of the cosmological constant, 
Qa > 0.7, lepton asymmetry must be small Jl5|— [l7[] , and 
T7 rsO.14. 



The vv annihilation mean free path for energies 
above the Z pole is much greater than the Hubble dis- 
tance H~ l . Therefore, those neutrinos whose present 
energy is E xes ± SE travel unabsorbed to red shifts 
z a = 2T Z /M Z = 0.06, which corresponds to distance 
280 Mpc(0.65//i). Within this distance the mean free 
path is A = l/(er anniZ 7V rc ii C ) = 5.3r7 _1 x 10 28 cm, where 
fann.z = 4nG F /y/2. Thus, within 50 Mpc from the ob- 
server a fraction 

p = 3 - 8 x 10 " (oil) (3) 

of the ultrahigh-energy neutrinos annihilate and produce 
"Z-bursts" . The nucleons and photons from Z decays 
undergo a cascade of scatterings off the background pho- 
tons but their energies remain largely above the GZK 
cutoff m . The resulting spectrum of these UHECR will 
be discussed in detail in an upcoming paper ]lq ]. Tak- 
ing into account the attenuation and absorption, each 
annihilation event produces N w 10 photons and pro- 
tons ||] whose energies exceed the GZK cutoff. The flux 
of UHECR beyond the GZK cutoff is predicted to be 




The mass of the X particle must exceed twice the res- 
onance energy E ICS = M z /2m v = 5.9(0.07eV/m„) x 
10 13 GeV. According to Rcf. @], in the mass range 
m x > 10 12 GeV, heavy particles decaying into neutri- 
nos are allowed by the present data as long as 

— P.' 1 > 2.4 x 10 5 

As pointed out by G. Sigl (private communication), 
the limit set by the Energetic Gamma Ray Experimncnt 
Telescope (EGRET) on the diffuse low-energy gamma- 
ray flux resulting from the Z-bursts could constrain 
this model. This question deserves further study. The 
other bounds usually mentioned in connection with the 
ultrahigh-energy neutrinos H do not apply because the 
spectrum of neutrinos in equation (Q) is peaked near m x . 
This differs drastically from the inverse power-law spec- 
trum assumed in Refs. M for astrophysical sources. In 
our case, the number of low-energy neutrinos is negligi- 
ble, and the only upper bound on the neutrino flux comes 
from the observations of atmospheric showers fl4[ ] . 

The flux of super-GZK cosmic rays in equation (|J) is 
consistent with the data for some range of parameters 
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around r x ~ (10 5 - 10 7 )t , m x ~ (10 13 - 10 15 ) GeV, 
m v ~ (0.01 — 1.0) eV. It is particularly interesting that 
the UHECR puzzle can be explained without fine-tuning 
the masses and the lifetimes of the hypothetical X par- 
ticles. The mass m x has to be at the scale where one 
generally expects to find some new physics. The same 
weakly interacting X particles could be the cold dark 
matter p2[ . Another interesting feature of our scenario 
is that the neutrino mass can be in the SK range. If the 
neutrino density is enhanced by either gravitational clus- 
tering Q (for neutrino masses ~ eV) or due to a lepton 
asymmetry of the universe ||, the parameter space for 
m x , t x increases accordingly. 

Since there is no natural reason for m x to be exactly 
equal to 2i? res , the neutrinos produced by the X particle 
decays in the local cluster of galaxies give a negligible 
contribution to the measured flux of cosmic rays. How- 
ever, those particles that decayed at red shift 



m x _ , = m x m v _ 
2E Ies M% 



(G) 



have just the right energy at present for an efficient anni- 
hilation. The arrival directions of UHECR point, there- 
fore, to their sources located in a thin spherical shell 
at red shift z e ± Sz, where 8z = (1 + z e ){T z /A4 z ) = 
0.03(1 + z e ). 

If m x ~ 10 14 GeV for m u — m SK , then z e ~ 0.1 — 1. In 
this case, one may hope to identify the sources with the 
remote clusters of galaxies at red shift z e . In principle, 
this would allow one to measure m x . The present data 
already shows some directional clustering of events [|l0| . 
There is a 1% probability for these correlations to be 
accidental jl0| . In the future, their statistical significance 
can be tested on a larger data sample. Recent studies of 
anisotropics in the UHECR jll| concentrated on the dark 
matter distribution in our galaxy and the local group of 
galaxies. In view of our results, a similar analysis of 
possible correlations with the known clusters of galaxies 
at some fixed (but yet unknown) red shift could reveal 
the locations of distant sources and help determine m x . 

In any case, our scenario does not predict any corre- 
lation of the directions of UHECR with the distribution 
of local dark matter. If the X particle decays are many- 
body, Sz is increased because of a wider spectrum of the 
emitted high-energy neutrinos. 

Non-thermal production of superheavy relic particles 
at the end of inflation J§|;[l2| and their role in cosmol- 
ogy [^9) have been the subject of intense studies re- 
cently. It is clear that £l x is a model-dependent param- 
eter, which can be ~ 0.1 in many realistic models H. 

We have presented a plausible candidate source of 
ultrahigh-energy neutrinos that can be the explanation 
of the cosmic rays beyond the GZK cutoff. Weakly in- 
teracting massive particle with mass m x > 10 13 GeV 
and lifetime t x > 10 5 io can decay into very energetic 
neutrinos. The expansion of the universe red shifts the 



energies of the neutrinos produced in a spherical shell at 
z = (m x m u /M z — 1) to the value for which the cross 
section of vv annihilation near Earth is large. These 
annihilations can create a flux of photons and nucleons 
with energies above the GZK cutoff that is sufficient to 
explain the present data. The possible directional corre- 
lations of UHECR events produced in this manner reflect 
the matter distribution at red shift z e ±5z. Our scenario 
is consistent with the neutrino masses in the range re- 
ported by the Super-Kamiokande experiment. 

This work was supported in part by the US Depart- 
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